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Comparison of Sc 3 and Y 3 
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Abstract 

CASSCF /CCI calculations are presented for the low-lying states of Y 3 . Compar- 
ison of the wave functions for Y 3 and SC 3 indicates substantial 4d-5p hybridization 
in Y 3 , but little 3d-4p hybridization in SC 3 . The increased 4d-5p hybridization 
leads to stabilization of 4 d 7 r bonding with respect to 4da bonding for equilateral 
triangle Y 3 , and also leads to 4d-5p bonding for linear geometries. These effects 
lead to a different ordering of states for equilateral triangle geometries and a smaller 
excitation energy to the linear configuration for Y 3 as compared to SC 3 . 
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I. Introduction 

Recently Knight, Woodward, VanZee, and Weltner[l] studied the ESR spectra 
of Sc 3 and Y 3 in rare gas matrices. The ESR spectrum for Sc 3 shows equivalent 
Sc atoms, which is consistent with an equilateral triangle ground state geometry. 
Calculations by Walch and Bauschlicher[2] find a 2 Aff ground state for Sc 3 with the 
configuration: 

4sa'j 2 4se' 4 3 da' 2 Sda" 1 

Here the 3da'j and 3da" orbitals are 3d bonding orbitals derived from the atomic 
3d<7 and 3d7r" orbitals, respectively. (The symmetry designations here are for D 3 j» 
symmetry and correspond to a a axis along a line connecting a given atom with the 
center of the molecule.) From a comparison of the Sc 3 and Ca 3 potential curves, it 
is apparent that these 3d bonding orbitals contribute significantly to the bonding 
in this molecule. Also, some 4p involvement due to 4s-4p hybridization is found in 
the bonding in both Ca 3 and Sc 3 . However, there is little admixture of 4p into the 
3d like bonding orbitals. 

For Y 3 , the ESR spectrum shows only two equivalent Y atoms, which is consistent 
with a non - equilateral triangle structure (such as an obtuse angle bent molecule 
or a linear molecule). It was of interest to see if the calculations would show 
significant differences between the bonding in Sc 3 and Y 3 , which could account for 
these differences which are observed experimentally. We find the major difference 
between these molecules is the increased importance of 5p-4d hybridization in the 
bonding orbitals of Y 3 , relative to 4p-3d hybridization in Sc 3 . This effect results 
in a different ordering of states for Y 3 as compared to Sc 3 and 5p-4d hybridization 
effects lead to a low-lying linear configuration for Y 3 , while the calculations on Sc 3 
clearly indicate an equilateral triangle ground state. In section II we briefly describe 
the calculational details. Section III describes the results obtained, and Section IV 
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presents the conclusions. 

II. Calculational Details 

The calculations are similar to the calculations on SC3 which have been reported 
previously. The calculations use a relativistic effective core potential (ECP) devel- 
oped by Hay and Wadt[3] which replaces the A r core, but includes the 4s and 4p 
core levels in the calculation along with the 5s, 5p, and 4d valence electrons. The Y 
basis set is the same as the basis set reported by Hay and Wadt, except that the 5p 
functions are replaced by the functions recommended by Walch, Bauschlicher, and 
Nelin[4] multiplied by 1.5 to make them suitable for describing 5s — >5p correlation. 
The s and p basis sets are contracted (2111) based on the atomic 4s and 4p orbitals, 
respectively. The d basis set was contracted (211) and a single set of 4f functions 
was added as a two-term fit to a Slater orbital(f = 1.40) [4]. The basis set and ECP 
were tested by comparing ECP and all-electron calculations [7] for the 5 £~ state 
of the Y2 molecule. The differences in the calculated potential curve parameters 
(ECP minus all-electron) were R e = -0.01 a 0 , w e = +4 cm -1 , and D e = +0.04 eV. 

The complete active space self consistent field (CASSCF) calculations^] included 
the 5s-5p and 4d derived orbitals in the active space. For near equilateral triangle 
geometries the active orbitals were a'j and e' from 5s, and e' from 4dcr, and a 2 
and e" from 4d7r". For linear geometries the active space included three o g , two 
o u , and one each n xu , n xg , 7r yu , and n yg . The nature of these orbitals is discussed 
in section III. As for the SC3 calculations^], the number of electrons in symme- 
try subgroups was constrained to be the same as in the dominant configuration. 
For C 2v geometries, the symmetry subgroups were i) aj and b2 and ii) a2 and bi, 
i.e. the orbitals were partitioned according to whether they were symmetric or 
asymmetric with respect to reflection in the molecular plane. For linear geome- 
tries, the symmetry subgroups were a, 7r x , and 7r y . The CASSCF calculations were 


3 


followed by contracted configuration interaction (CCI) calculations^]. Here the ref- 
erence space included those configurations with coefficients greater than 0.1 in the 
CASSCF wavefunction near R e (five, ten, and thirteen reference configurations for 
2 Aj, 2 E', and 2 E+, respectively). This choice of reference configurations leads to 
valence Cl energies within 0.1 eV of the CASSCF energy and reference percentages 
in the wavefunction of about 90 % . The calculations were carried out with the 
MOLECULE[5]- SWEDEN[6] system of programs. 

III. Discussion 

Table I summarizes the dominant configurations, geometries, symmetric stretch 
frequencies, T e ’s, and Mulliken populations for the low-lying states of SC3 and Y3. 
Tables II and III give CASSCF and CCI energies for equilateral triangle geometries 
for the 2 Aj, 2 A2, and 2 E' states of Y3. 

The SC3 molecule has a 2 A2 ground state with the configuration: 

4sa'j 2 4se' 4 3da'j 2 Zda,'^ 1 

and a very low-lying 2 k\ state arising from the same limit with the configuration: 
4sai 2 4se' 4 Sda'j 1 3da2 2 

Formally these configurations arise from three Sc atoms in the 4s 2 3d 1 ground state. 
There is also a low-lying 2 E' excited state arising from two ground state Sc atoms 
and one Sc atom in the 4s 1 3d 2 excited state: 

4sa'j 2 4se' 3 3da( 2 3da' 2 ' 2 

From Table I there are several obvious differences between Y3 and SC3. The first 
difference is the larger p population for Y3 as compared to SC3. Here Y3 has a 5p 
population of 2.25 for the 2 Aj ground state, while SC3 has a 4p population of 1.21 for 
the 2 Aj state. For SC3 it was argued that the 4p population involved admixture of 
terms derived from the 4s 1 4p 1 3d 1 Sc atomic configuration, which is at an excitation 
energy of 1.96 eV[lO]. For Y atom the Ss^pMd 1 configuration is at 1.91 eV.10]; 
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thus, one would expect similar p populations if this were the only mechanism for 
mixing in p character. However, for the Y atom the lowest excited state has the 
configuration 5s 2 5p J and is at an excitation energy of 1.33 eV[lO]. By comparison, 
the analagous state for Sc atom is not observed experimentally, but is calculated to 
be at 3.06 eV from SCF calculations. Given the low-lying 5s 2 5p J state of Y atom 
one expects that the ground state of Y 3 could arise from two atoms in the 5s 2 4d 1 
ground state and one atom in the low-lying 5s 2 5p 1 state. States derived from this 
limit would still mix in 5s 1 5p 1 4d 1 character leading to p populations about 1.0 
greater than for states arising from three ground state atoms. 

A second difference between SC 3 and Y 3 is the different ordering of states. For SC 3 
the 2 Aj , 2 A 2 , and 2 E' states are all close in energy; while, for Y 3 the ground state is 
2 A' and the 2 Aj — + 2 Aj excitation energy is 0.75 eV (from CASSCF calculations) 
for Y 3 as compared to -0.07 eV for SC 3 . This implies a stabilization of the 4da 2 
orbital relative to the 4da'j orbital for Y 3 . The reason for this is obvious from the 
populations of the 4da£ and 4dai orbitals. From Table IV. and Fig. 1 one sees 
a considerable admixture of 5p-4d character in the 4d bonding orbitals of Y 3 , but 
much less admixture of 4p-3d character for SC 3 . This difference arises because of 
the more comparable 5p-4d orbital sizes as compared to the 4p-3d orbitals[7] and 
the smaller excitation energy to the 5s 2 5p 1 state of Y atom as compared to the 
corresponding state of the Sc atom. Nonetheless, the 5p orbitals are still larger 
than the 4d orbitals so that distances corresponding to good 4d-4d overlaps [7] are 
well inside the optimal bonding radius for the 5p orbitals. Furthermore, the overlap 
integral Ss pir -$ pn increases monotonically toward one as R goes toward zero, but the 
overlap integral Sspe- 5 P a decreases inside the optimal bonding radius, eventually 
goes through zero, and becomes minus one for R equal zero. Thus, 4d-5p admixture 
is more favorable for 5 p 7 r than for 5p o and this effect appears to be responsible for 
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the stabilization of the 4da" orbital relative to the 4da x orbital for Y3. 

We now consider non-equilateral triangle geometries. For C2 V symmetry the 2 A' X , 
2 A' 2 ', and 2 E' states correspond to 2 A X , 2 Bi, and 2 2 Ai and 2 B2, respectively. Here 
the configurations are: 

2 Aj = la 2 2a 2 lb 2 lb 2 3a} 

2 Bj = la x 2a x lb 2 lb x 3a x 
2 2 Ai = la^aJlbjlbjSaf 
2 B 2 = la 2 2a?lb*lb 2 3a 2 

We expect the 2 A x and 2 A2 states to have equilateral triangle geometries, while 
the 2 E' state will Jahn- Teller distort leading to an acute angle minimum for the 
2 B2 component and an obtuse angle minimum for the 2 Aj component. These 
expectations are confirmed by calculations. Table V gives CASSCF energies for the 
2 Ai state of Y3 for non-equilateral triangle geometries (with ry_y fixed at 6.3 ao 
which is the optimal value for the 2 A X state from Table II.). Here one sees that the 
energy of the 2 Ai state goes up for obtuse angle geometries. Table V also shows 
energies for the 2 2 Ai state obtained as a second root of the CASSCF. Here one sees 
that, as expected, this state favors obtuse angle geometries. 

For linear geometries, the 2 2 Aj configuration corresponds to: 

2 £+ = 10g2Cg Ia2l7r 2 u l7r 2 u 

Figure 2 shows key bonding orbitals for Y3 for linear geometries. The 2 o g orbital is 
basically a localized 4d a orbital. The ln g orbital is a bonding orbital of primarily 
4 dn character, while the l?r u orbital is of mixed character( mostly 5p7r on the center 
atom and 4dtf on the end atom). The occupations of these orbitals for the 2 E+ 
state are 2 a g = 1.01, l7r u = 2 x 1.11 and l7r ff = 2 x 0.89. Given the very similiar 
populations for the l7r u and l7r p orbitals, this state is best thought of as an open 
shell doublet state derived from: 
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6 £+ = lCg20gl0^l7rl u l*lg\*\ U \*\g 

This change in the character of the wavefunction may be thought of in terms of a 
curve crossing between the 2 2 Ai configuration which is stable for near equilateral 
triangle geometries and the 2 £+ state which is stable for linear geometries. However, 
the 2 £+ state does retain some of the character of the 2 E' state as evidenced by 
the 4d population of 3.82. 

Table VI gives CASSCF and CCI energies for the 2 £+ state of Y3 for linear 
geometries. Table VII compares Mulliken populations for the 2 o g , l7r IU , and l 7 r xg 
orbitals of the 6 £+ states of Y3 and SC3, while plots of these orbitals are shown in 
Fig. 2. (The 6 £+ orbitals are found to be very similar to the 2 £+ orbitals for Y3.) 
From Table VII and Fig. 2, one sees much more admixture of 5p and 4d in the 7r 
bonding orbitals of Y3 as compared to the corresponding orbitals of SC3. This is 
particularly evident for the l7r xu orbital. From Table VII it is seen that the latter 
orbital has about equal parts of 5p and 4d character for Y3, but is mostly 4p like for 
SC3. These differences are also evident in the plots of these orbitals in Fig. 2. The 
l7r IU orbital of Y3 shows more 4d contours which is indicative of a Y(5p7r)-Y(4d7r) 
bond. The increased 5p-4d bonding in Y 3 is reflected in the small excitation energy 
to the 6 £+ state of 0.31 eV for Y3 as compared to an excitation energy of 1.28 eV 
for the corresponding excitation in Scs( see Table I.). 

The ESR spectrum of Y3 is consistent with non equivalent Y atoms. The lowest 
state in the calculations is 2 Aj, which is expected to have an equilateral triangle 
geometry and is not consistent with the experiment (see Table V). It would be 
tempting to argue that the ground state of Y3 is the linear 2 £+ state. This state 
is fully consistent with the ESR spectrum, since only parallel and perpendicular 
components of the g and A tensors were needed to simulate the spectrum. Also 
the singly occupied orbital has only about 5 % 5s character. The small amount of 
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5s character is consistent with the primarily 5p and 4d character of the 2 cr p , l 7 r xu , 
l7r yu , and l7T vfl orbitals which are coupled together into a doublet spin state. 
The calculations do show that linear configurations are lower for Y 3 than for Sc 3 , 
but the computed excitation energy of 1.54 eV is quite large. 

It is not clear whether the computed excitation energy to the 2 E+ state could 
be off by more than 1.5 eV. We do note that for the ground state a E+ is 
computed to be 0.40 eV higher than the open shell 7 E+ state. Similarly, for Y 2 , 
computations show a 5 E~ ground state, but ESR does not observe such a state and 
one possibility is that the true ground state is which is computed to be 0.87 eV 
higher. Thus, there is reason to believe that the calculations could have significant 
errors in the excitation energies for states of very different character. For the cases 
above we expect that the calculation will be biased toward open shell high-spin 
states over closed shell states. Similarly, for Y 3 the linear 2 E^ state has a much 
less compact zero order wavefunction than do the equilateral triangle states which 
are SCF dominated. Thus, one might also expect significant errors in the excitation 
energies here, and the ground state of Y 3 may be linear in agreement with one 
interpretation of the ESR spectrum of Y 3 . Finally, we note that although the CCI 
calculations for the 2 E+ state correspond to more than one million configurations, 
the reference space is only about 90 % of the Cl wavefunction. This problem, which 
has proven to be typical of calculations on transition metal dimers and trimers, 
is thought to arise because of the approximation of not including explicit nd shell 
correlation in the CASSCF wavefunction. 

In spite of the difficulties in computing the excitation energies for these systems, 
these calculations show significant qualitative differences in the bonding in these 
molecules. This qualitative understanding of the bonding is important, since very 
little is known about the electronic structure of small transition metal clusters. Also 
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based on experience for transition metal diatomics [7], we expect that the computed 
geometries and force constants will be accurate and may provide some guidance to 
experimental studies on these molecules. 

IV. Conclusions 

We conclude that the bonding in Sc 3 and Y 3 differ mainly in terms of a larger p 
involvement in the bonding in Y 3 as compared to SC 3 . This larger p involvement 
arises because the 5s 2 5p ! state of the Y atom is low-lying(the first excited state), 
while the corresponding state for the Sc atom is at a large excitation energy. The 
result is that the lowest state of Y 3 for equilateral triangle geometries arises from 
two Y atoms in the 5s 2 4d* state and one Y atom in the 5s 2 5p 1 state, leading to a 
5p population greater than two. For Sc 3 , on the other hand, the lowest state arises 
from three ground state atoms and has a 4p population of about one. 

In addition to the larger p population for Y 3 , the 5p and 4d orbitals are more 
comparable in size than the 4p and 3d orbitals, which leads to significant 5p-4d 
hybridization in Y 3 , but 4p-3d hybridization is not important in Sc 3 . The 5p-4d 
hybridization stabilizes the orbital relative to the 4 da'j( 4 da) orbital. 

The latter effect leads to a 2 A\ ground state for equilateral triangle Y 3 , wheras, 
equilateral triangle SC 3 has a 2 A 2 ground state. 

For linear geometries Y 3 has a low-lying 2 £+ state. This state is much lower for 
Y 3 than for Sc 3 because of the formation of 5p-4d n bonds in Y 3 . If this state were 
the ground state of Y 3 it would be fully consistent with the matrix ESR spectrum 
of Y 3 . This assignment would require that the computed excitation energy to this 
state be in error by more than 1.5 eV. However, large errors in excitation energies 
for transition metal systems for states of different character have been observed for 
other systems, and thus the ground state of Y 3 may be linear in agreement with 
one possible interpretation of the ESR spectrum. 
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Table I. Configurations, Potential Surface Parameters, and Mulliken Populations 
for SC3 “ and Y3. 
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Table II. CASSCF Energies for Y 3 (Equilateral Triangle Geometries) 


R(Y-Y) 

2 a; 

2 E / 

2 A£ 

5.25 


-112.37860 


5.50 

-112.40823 

-112.40455 


5.75 


-112.41844 


6.00 

-112.44621 

-112.42295 

-112.41858 

6.30 

-112.45092 

-112.41944 


6.50 

-112.44918 



7.00 

-112.43454 
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Table III. CCI Energies for Y? (Equilateral Triangle Geometries) 


R(Y-Y) 

2 A'i 

2 E' 

5.50 


-112.56399 

5.75 


-112.57800 

6.00 

-112.5 >0 

-112.58221 

6.30 

-112.59c 97 

-112.57768 

6.50 

-112.59524 


7.00 

-112.57640 
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Table IV. Mulliken Populations for Selected Orbitals for SC 3 and Y 3 for the 2 A\ 
state. 



s 

P 

d 

f 

Sc 3 (3dai) 

0.01679 

0.07461 

0.73747 

0.01106 

Y s (4da') 

0.03922 

'.23463 

0.71649 

0.02088 

Sc 3 (3da' 2 ') 

0.00000 

0.18008 

1.39599 

0.01473 

Y 3 (4da 2 ) 

0.00000 

0.48791 

1.34106 

0.01985 


Table V. CASSCF Energies for Y 3 (Non-Equilateral Triangle Geometries) 


R(Y-Y) 

e 

2 Ai 

2 2 Ai 

6.30 

60. 

-112.45092 

-112.39422 

6.30 

75. 

-112.43708 

- 112.40411 

6.30 

90. 

-112.41910 

-112.39700 

6.30 

120 . 

-112.3985 

-112.3734 
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Table VI. CASSCF and CCI Energies for Y 3 (Linear Geometries) 


R(Y-Y) 

CASSCF 

CCI 

5.50 

-112.41611 

-112.53973 

5.75 

-112.42374 

-112.54227 

6.00 

-112.42587 

-112.54239 

6.30 

-112.42279 

-112.54006 
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Table VII. Mulliken Populations for Selected Orbitals for SC 3 and Y 3 for the 6 E + 


state. 


Sc 3 (2oj,) 

s 

0.03671 

P 

0.00218 

d 

0.98366 

f 

0.00001 

Y 3 (2o 9 ) 

0.25686 

0.06009 

0.67001 

0.00437 

Sc 3 (l7riu) 

0.00000 

0.81349 

0.17874 

0.00868 

Y 3 (l7riu) 

0.00000 

0.45445 

0.52200 

0.02225 

Sc 3 (l7r* p ) 

0.00000 

0.10594 

0.88688 

0.00631 

Y 3 (l7T r9 ) 

0.00000 

0.13547 

0.85645 

0.00938 


18 



Figure Captions. 


Figure 1. Contour plots of selected orbitals of the 2 A 2 state of Sc 3 and the 2 Aj 
state of Y 3 for equilateral triangle geometries. The orbitals are plotted in a plane 
perpendicular to the molecular plane and passing through one atom and the center 
of the opposite bond. Positive contours are denoted by solid lines, negative contours 
are denoted by dotted lines, and nodal surfaces are denoted by dashed lines. 

Figure 2. Contour plots of selected orbitals for the 6 £+ state of Sc 3 and Y 3 for 
linear geometries. The orbitals are plotted in the molecular plane. The conventions 
are the same as in Fig. 1. 










